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Abstract 

I work out the Newtonian and general relativistic effects due to an 
isotropic mass loss M /M of a body on the orbital motion of a test par- 
ticle around it; the present analysis is also valid for a variation G/G 
of the Newtonian constant of gravitation. Concerning the Newtonian 
case, I use the Gauss equations for the variation of the elements and 
obtain negative secular rates for the osculating semimajor axis a, the 
eccentricity e and the mean anomaly Ai, while the argument of peri- 
center w does not experience secular precession; the longitude of the 
ascending node 17 and the inclination i remain unchanged as well. The 
anomalistic period is different from the Keplerian one: it turns out 
to be larger than it. The true orbit, instead, expands, as shown by a 
numerical integration of the equations of motion with MATHEMAT- 
ICA; in fact, this is in agreement with the seemingly counter-intuitive 
decreasing of a and e because they refer to the osculating Keplerian 
ellipses which approximate the trajectory at each instant. A com- 
parison with the results obtained with different approaches by other 
researchers is made. General relativity induces positive secular rates 
of the semimajor axis and the eccentricity completely negligible in the 
present and future evolution of the solar system. 
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1 Introduction 

In this paper I investigate the classical orbital effects induced by an isotropic 
variation M /M of the mass of a central body on the motion of a test par- 
ticle; my analysis is valid also for a change G/G oi the Newtonian constant 
of gravitation. This problem, although interesting in itself, is not only an 
academic one because of the relevance that it may have on the ultimate 
destiny of planetary companions in many stellar systems in which the host 
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star experiences a mass loss, like our Sun [IJ. With respect to this aspect, 
my analysis may be helpful in driving future researches towards the imple- 
mentation of long-term N-body simulations including the temporal change 
of GM as well, especially over timescales covering paleoclimate changes, up 
to the Red Giant Branch (RGB) phase in which some of the inner planets 
should be engulfed by the expanding Sun. Another problem, linked to the 
one investigated here, which has recently received attention is the observa- 
tionally determined secular variation of the Astronomical Unit [21 [3l HI [5] . 
Moreover, increasing accuracy in astrometry pointing towards microarcsec- 
ond level p], and long-term stability in clocks [7] require to consider the 
possibility that smaller and subtler perturbations will be soon detectable in 
the solar system. Also future planetary ephemerides should take into ac- 
count M/M. Other phenomena which may, in principle, show connections 
with the problem treated here are the secular decrease of the semimajor axes 
of the LAGEOS satellites, amounting to 1.1 mm d~^, [8j and the increase 
of the lunar orbit's eccentricity of 0.9 x 10~^^ yr~^ [9]. 

Many treatments of the mass loss-driven orbital dynamics in the frame- 
work of the Newtonian mechanics, based on different approaches and laws 
of variation of the central body's mass, can be found in literature; see, e.g., 
[iniiniliailSllIllIiailSlIITliniElllland references therein. However, they 
are sometimes rather confused and involved, giving unclear results concern- 
ing the behavior of the Keplerian orbital elements and the true orbit. 

The plan of the paper is as follows. Section [2] is devoted to a theoreti- 
cal description of the phenomenon in a two-body scenario. By working in 
the Newtonian framework, I will analytically work out the changes after 
one orbital revolution experienced by all the Keplerian orbital elements of a 
test particle moving in the gravitational field of a central mass experiencing 
a variation of its GM linear in time. Then, I will clarify the meaning of 
the results obtained by performing a numerical integration of the equations 
of motion in order to visualize the true trajectory followed by the planet. 
Concerning the method adopted, I will use the Gauss perturbation equa- 
tions \19\ [20] , which are valid for generic disturbing accelerations depending 
on position, velocity and time, the "standard" Keplerian orbital elements 
(the Type I according to, e.g., Ref. [1^) with the eccentric anomaly E as 
"fast" angular variable. Other approaches and angular variables like, e.g. 
the Lagrange perturbation equations [l9l[20], the Type II orbital elements 
[16j and the mean anomaly M could be used, but, in my opinion, at a price 
of major conceptual and computational difficultie^. With respect to possi- 

^ Think, e.g., about the cumbersome expansions in terms of the mean anomaly and the 
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ble connections with realistic situations, it should be noted that, after all, 
the Type I orbital elements are usually determined or improved in standard 
data reduction analyses of the motion of planets and (natural and artificial) 
satellites. Instead, my approach should, hopefully, appear more transpar- 
ent and easy to interpret, although, at first sight, some counter-intuitive 
results concerning the semimajor axis and the eccentricity will be obtained; 
moreover, for the chosen time variation of the mass of the primary, no ap- 
proximations are used in the calculations which are quite straightforward. 
However, it is important to stress that such allegedly puzzling features are 
only seemingly paradoxical because they will turn out to be in agreement 
with numerical integrations of the equations of motion, as explicitly shown 
by the numerous pictures depicted. Anyway, the interested reader is advised 
to look also at Ref. [16) for a different approach. In Section [3] I will work 
within the general relativistic gravitoelectromagnetic framework by calcu- 
lating the gravitoelectric effects on all the Keplerian orbital elements of a 
freely falling test particle in a non-stationary gravitational field. Section U] 
is devoted to a discussion of the findings of other researchers and contains 
some numerical calculations concerning the previously mentioned orbital 
phenomena of LAGEOS and the Moon. Section [5] summarizes my results. 



2 Analytical calculation of the orbital effects by 

By defining 

11 = GM (1) 

at a given epoch to, the acceleration of a test particle orbiting a central body 
experiencing a variation of ^ is, to first order in t — Iq, 



ii 



1 + 



{t - to) 



2) 



with fi = fj-\t=to- A "will be assumed constant throughout the temporal 
interval of interest At = t — to, as it is, e.g., the case for most of the 
remaining lifetime of the Sun as a Main Sequence (MS) star Note that /i 
can, in principle, be due to a variation of both the Newtonian gravitational 
constant G and the mass M of the central body, so that 

(3) 



M 

G^M' 



Hansen coefficients, the subtleties concerning the choice of the independent variable in the 
Lagrange equations for the semimajor axis and the eccentricity |19j . 
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Moreover, while the orbital angular momentum is conserved, this does not 
happen for the energy. 

By limiting ourselves to realistic astronomical scenarios like our solar 
system, it is quite realistic to assume that 



(t - to) « 1 (4) 



over most of its remaining lifetime: indeed, since M/M is of the order o:f 
10~^^ yr"^ for the Sun [T], the condition ([3]) is satisfied for the remaining^ 
~ 7.58 Gyr before the Sun will approach the RGB tip in the Hertzsprung- 
Russell Diagram (HRD). Thus, I can treat it perturbatively with the stan- 
dard methods of celestial mechanics. 

The unperturbed Keplerian ellipse at epoch to, assumed coinciding with 
the time of the passage at perihelion tp, is characterized by 

r = a(l — e cos E), 
dt= [^-=^f^)dE, 



cos/ 
sin / : 



cos E—e 



1—e cos E ' 

yjl—e^ sing 
1— ecosi? ' 



where a and e are the semimajor axis and the eccentricity, respectively, 
which fix the size and the shape of the unchanging Keplerian orbit, n = 
is its unperturbed Keplerian mean motion, / is the true anomaly, 
reckoned from the pericentre, and E is the eccentric anomaly. This would 
be the path followed by the particle for any t > tp\i the mass loss would 
suddenly cease at ty. Instead, the true path will be different because of 
the perturbation induced by /i and the orbital parameters of the osculating 
ellipses approximating the real trajectory at each instant of time will slowly 
change in time. 

^About 80% of such a mass-loss is due to the core nuclear burning, while the remaining 
20% is due to average solar wind. 

^The age of the present-day MS Sun is 4.58 Gyr, counted from its zero-age MS star 
model [T]. 



4 



2.1 The semimajor axis and the eccentricity 

The Gauss equation for the variation of the semimajor axis a is |19L [20] 



da 



n 



eAr sin f + Ar 



(6) 



where Ar and Ar are the radial and transverse, i.e. orthogonal to the di- 
rection of f, components, respectively, of the disturbing acceleration, and 
p = a{l — e^) is the semilatus rectum. In the present case 



(t-tp), 



(7) 



i.e. there is an entirely radial perturbing acceleration. For /i < 0, i.e. a 
decrease in the body's GM, the total gravitational attraction felt by the 
test particle, given by ([2]), is reduced with respect to the epoch tp. In order 
to have the rate of the semimajor axis averaged over one (Keplerian) orbital 
revolution ([7|) must be inserted into ([6]), evaluated onto the unperturbed 
Keplerian ellipse with ([5]) and finally integrated over ndt/2TT from to 27r 
because n/27r = (see below). Note that, from ([5]), it can be obtained 



t-t. 



E — e sin E 



n 



As a result, I have 



da 

la 



e 
vr 



27r 



{E — e sin E) sin E 
(1 — e cos E)"^ 



dE 



e 

1-e 



(8) 



(9) 



Note that if ^ decreases a gets reduced as well: (a) < 0. This may be 
seemingly bizarre and counter-intuitive, but, as it will be shown later, it is 
not in contrast with the true orbital motion. 

The Gauss equation for the variation of the eccentricity is p!9| 



de 
~dt 



na 



For A = Ar, it reduces to 



so that 



Ar sin / + Ar 



de 
dt 



de 
di 



2ae 



(1 + e) 



cos / + - ( 1 



da 
A 



(10) 



(11) 



(12) 
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also the eccentricity gets smaller for /t < 0. 

As a consequence of the found variations of the osculating semimajor 
axis and the eccentricity, the osculating orbital angular momentum per unit 
mass, defined by = fj,a{l — e^), remains constant: indeed, by using Q 
and (fT2]) . it turns out 

'dL^ \ 

fi {d) {I - e"^) - 2fiae {e) = 0. (13) 



dt 

The osculating total energy £ = —fi/2a decreases according to 

^(d) = f-^U. (14) 



dt I Ic? \ 1 — e / a 

Moreover, the osculating Keplerian period 



n 




pKep^'^^27rJ-, (15) 



which, by definition, yields the time elapsed between two consecutive per- 
ihelion crossings in absence of perturbation, i.e. it is the time required to 
describe a fixed osculating Keplerian ellipse, decreases according to 

/dP^\ ^ 3 KepM ^ 6^eA /ay/2 

\ dt / 2 a {l-e)\fij ' ^ ^ 

As I will show, also such a result is not in contrast with the genuine orbital 
evolution. 



2.2 The pericentre, the node and the inchnation 

The Gauss equation for the variation of the pericentre cj is \19\ [20] 

sin/ 



du) 
~dt 



nae 



-Ar cos / + 1 + 



P 



.dn 



cos I 



dt 



(17) 



where i and Q are the the inclination and the longitude of the ascending 
node, respectively, which fix the orientation of the osculating ellipse in the 
inertial space. Since dQ/dt and di/dt depend on the normal component 
Ai, of the disturbing acceleration, which is absent in the present case, and 
A = Ar, I have 

"2^ (^ . 



duj 
~di 



27re 



esin£')(cos£' — e) 



dE = 0: 



(18) 



(1 — e cos E)"^ 

the osculating ellipse does not change its orientation in the orbital plane, 
which, incidentally, remains fixed in the inertial space because Ai, = and, 
thus, dVL/dt = di/dt = 0. 
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2.3 The mean anomaly 



The Gauss equation for the mean anomaly Ai, defined as = n{t — tp), 
[HIISO] is 

dM 2 r r ^/dw .d^\ 

——=n Ar Vl-e^ — + COS2— . (19) 

di na a \ dt dt j 

It turns out that, since 

2 r 2/y 

Ar-dt = (E-e sin E)dE, (20) 

na a n-^a-^ 

the mean anomaly changes uniformly in time at a slower rate with respect 
to the unperturbed Keplerian case for /i < 0. 



2.4 Numerical integration of the equations of motion and 
explanation of the seeming contradiction with the ana- 
lytical results 

At first sight, the results obtained here may be rather confusing: if the 
gravitational attraction of the Sun reduces in time because of its mass loss 
the orbits of the planets should expand (see the trajectory plotted in Figure 
m numerically integrated with MATHEMATICA) , while I obtained that the 
semimajor axis and the eccentricity undergo secular decrements. Moreover, 
I found that the Keplerian period P^'^p decreases, while one would expect 
that the orbital period increases. 

In fact, there is no contradiction, and my analytical results do yield us 
realistic information on the true evolution of the planetary motion. Indeed, 
a, e and P^'^p refer to the osculating Keplerian ellipses which, at any instant, 
approximate the true trajectory; it, instead, is not an ellipse, not being 
bounded. Let us start at tp from the osculating pericentre of the Keplerian 
ellipse corresponding to chosen initial conditions: let us use a heliocentric 
frame with the x axis oriented along the osculating pericentre. After a 
true revolution, i.e. when the true radius vector of the planet has swept 
an angular interval of 27r, the planet finds itself again on the X axis, but 
at a larger distance from the starting point because of the orbit expansion 
induced by the Sun's mass loss. It is not difficult to understand that the 
osculating Keplerian ellipse approximating the trajectory at this perihelion 
passage is oriented as before because there is no variation of the (osculating) 
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y, AU 




Figure 1: Black continuous line: true trajectory obtained by numerically 
integrating with MATHEMATICA the perturbed equations of motion in 
Cartesian coordinates over 2 yr; the disturbing acceleration ([2]) has been 
adopted. The planet starts from the perihelion on the x axis. Just for 
illustrative purposes, a mass loss rate of the order of 10~^ yr"^ has been 
adopted for the Sun; for the planet initial conditions corresponding to a = 1 
AU, e = 0.8 have been chosen. Red dashed line: unperturbed Keplerian 
ellipse at t = Iq = tp. Blue dash-dotted line: osculating Keplerian ellipse 
after the first perihelion passage. As can be noted, its semimajor axis and 
eccentricity are clearly smaller than those of the initial unperturbed ellipse. 
Note also that after 2 yr the planet has not yet reached the perihelion as 
it would have done in absence of mass loss, i.e. the true orbital period is 
longer than the Keplerian one of the osculating red ellipse. 
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argument of pericentre, but has smaller semimajor axis and eccentricity. 
And so on, revolution after revolution, until the perturbation theory can be 
applied, i.e. until — tp) « 1. In Figured] the situation described so 

far is qualitatively illustrated. Just for illustrative purposes I enhanced the 
overall effect by assuming fi/ fj, ~ 10^^ yr^^ for the Sun; the initial conditions 
for the planet correspond to an unperturbed Keplerian ellipse with a = 1 AU, 
e = 0.8 with the present-day value of the Sun's mass in one of its foci. It is 
apparent that the initial osculating red dashed ellipse has larger a and e with 
respect to the second osculating blue dash-dotted ellipse. Note also that the 
true orbital period, intended as the time elapsed between two consecutive 
crossings of the perihelion, is larger than the unperturbed Keplerian one of 
the initial red dashed osculating ellipse, which would amount to 1 yr for the 
Earth: indeed, after 2 yr the planet has not yet reached the perihelion for 
its second passage. 

Now, if I compute the radial change Ar{E) in the osculating radius vector 
as a function of the eccentric anomaly E I can gain useful insights concern- 
ing how much the true path has expanded after two consecutive perihelion 
passages. Prom the Keplerian expression of the Sun-planet distance 



one gets the radial component of the orbital perturbation expressed in terms 
of the eccentric anomaly E 



it agrees with the results obtained by, e.g., Casotto in Ref. |21) . Since 



r = a(l — e cos E) 



(22) 



Ar{E) = {1 - e cos E) Aa - a cos E Ae + ae sin E AE; 



(23) 



Aa 




< 



Ae 




(24) 



AE 



AM+sinE Ae \ _ i 
l—ecosi? / n 



I) [AiE) + B{E) + CiE)] 



with 



r A{E) 



_E^+2e(cos E-1) 



1— e cos E 



< 




(25) 



[ C{E) 



(1— e^) sin£;(sin B— e cos E) 
(l-ecos£;)2 
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Ar{E) 



[V{E)+:FiE)], 



(26) 



-2(_Sin - COS iij H 1-ecosE (l-ecosE)^ 



it follows 

with 
V{E) 

F[E) 

(27) 

From (I26p and (1270 it turns out that for E > Q Ar{E) never vanishes; 
after one orbital revolution, i.e. after that an angular interval of 27r has 
been swept by the (osculating) radius vector, a net increase of the radial 
(osculating) distance occurs according tc0 



l-e^ 



1— e cos J 



j < COS £'(sin E — E cos E) + sin E 



l+e-(l+e) cos_B-_Esin_B 
1— e cos E 



Ar(27r) - Ar(0) = Ar(27r) 



27r / ii 

a - 

n \fi 



(1-e). 



(28) 



This analytical result is qualitatively confirmed by the differenc^ Ar{t) 
between the radial distances obtained from the solutions of two numerical 
integrations with MATHEMATICA of the equations of motion over 3 yr 
with and without the initial conditions are the same. For illustrative 
purposes I used a = 1 AU, e = 0.01, = —0.1 yr~^. The result is depicted 
in Figure El Note also that and ^I7\i tell us that the shift at the aphelion 
is 

Arfvr) = l( ^-t-^ 1 Ar(27r), (29) 



2 VI - 

in agreement with Figure [U where it is 4.5 times larger than the shift at the 
perihelion. 

Since Figure [D tells us that the orbital period gets larger than the Kep- 
lerian one, it means that the true orbit must somehow remain behind with 
respect to the Keplerian one. Thus, a negative perturbation At in the 
transverse direction must occur as well; see Figure [3l 

Let us now analytically compute it. According to Ref. [21], it can be 
used 

asini? / . „ , , . ^ , , , 

+ a^/l-e^ AE + r{Auj + Ancosi). (30) 

V 1 — 



At 



* According to ^ and (EZj, Ar(0) = 0. 

^Strictly speaking, Ar and the quantity plotted in Figure [2] are diflerent objects, but, 
as the following discussion will clarify, I can assume that, in practice, they are the same. 
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Ar (AU) 




Figure 2: Difference Ar{t) between the radial distances obtained from the 
solutions of two numerical integrations with MATHEMATICA of the equa- 
tions of motion over 3 yr with and without fi/fi; the initial conditions are 
the same. Just for illustrative purposes a mass loss rate of the order of 
— 10~^ yr~^ has been adopted for the Sun; for the planet initial conditions 
corresponding to a = 1 AU, e = 0.01 have been chosen. The cumulative 
increase of the Sun-planet distance induced by the mass loss is apparent. 



11 



\ True orbit 

Unperturbed ellipse '"-^ 



\ 

\ 




Line of the 
apsides 



Figure 3: Radial and transverse perturbations Ar and At of the Keplerian 
radius vector (in blue); the presence of the transverse perturbation At 
makes the real orbit (in red) lagging behind the Keplerian one. 
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By recalling that, in the present case, AQ = and using 

l + e- (l + e)cos^-£;sin^' 



Auj 



ne \iJ, 



1 — e cos E 



(31) 



it is possible to obtain from and (|HT]) 
a / /i 



At{E) 
with 



n \fij (1 — e cos E) 

n{E) 

X{E) 
J{E) 



[G{E) +n{E)+ Z{E) + J{E) + IC{E)] 



(32) 



sin E(E cos E — sinii^), 
^^""^"'^^^ [(1 + e)(cos£; - 1) + ^sin^] , 
^2 + 2e(cos^- 1), 

(1— e^){e cos B— sin E) 



(33) 



IC{E) 



sinE 



l—e cos E 



(l+e)(l-ecos E)-E sin E 
l—e cos E 



From ([32]) and (j33|) it turns out that for E > At{E) never vanishes; at 
the (osculating) time of perihelion passage 



AT(27r) - At(0) 



Air" I 



— a 
n V /i 



l + e 
1 - e 



< 0. 



(34) 



This means that when the Keplerian path has reached the perihelion, the 
perturbed orbit is still behind it. Such features are qualitatively confirmed 
by Figure m 

From a vectorial point of view, the radial and transverse perturbations 
to the Keplerian radius vector r yield a correction 



so that 

The length of A is 



A = Ar f + Ar f, 
rpeit = r + A. 

A(^) = VAr(^)2 + Ar(E)2; 



(35) 
(36) 

(37) 
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([28]) and ([32]) tell us that at perihelion it amounts to 



A(2vr) = Ari2n)J 1 + ^^'j^^,- (38) 



The angle ^ between A and r is given by 



at perihelion it is 



tane(2vr) = -2vr^^^i±^, (40) 

i.e. ^ is close to —90 deg; for the Earth it is —81.1 deg. Thus, the difference 
6 between the lengths of the perturbed radius vector rpert and the Keplerian 
one r at a given instant amounts to about 

(5«Acos^; (41) 

if fact, this is precisely the quantity determined over 3 yr by the numerical 
integration of Figure El At the perihelion I have 



5 = Ar(27r)^/l + 4^2_ll±^cose; (42) 



since for the Earth 



Wl + 47r2iii^ cose = 1.0037, (43) 

it holds 

6^Ar{27r). (44) 

This explains why Figure [2] gives us just Ar. 

Since the approximate calculations of other researchers often refer to 
circular orbits, and in view of the fact that when a Sun-like star evolves into a 
giant tidal interactions circulariz^ the orbit of a planet [2_2] , it is interesting 
to consider also such limiting case in which other nonsingular osculating 
orbital elements must be adopted. The eccentricity and the pericentre lose 
their meaning: thus, it is not surprising that (I12|) . although formally valid 
for e — )• 0, yields a meaningless result, i.e. the eccentricity would become 
negative. Instead, the semimajor axis is still valid and ([9]) predicts that 
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Y, AU 




Figure 4: Black continuous line: true trajectory obtained by numerically 
integrating with MATHEMATICA the perturbed equations of motion in 
Cartesian coordinates over 2 yr ; the disturbing acceleration ((2j) has been 
adopted. The planet starts from a point on the x axis. Just for illustrative 
purposes, a mass loss rate of the order of — 10~^ yr~^ has been adopted 
for the Sun; for the planet initial conditions corresponding to a = 1 AU, 
e = 0.0 have been chosen. Red dashed line: unperturbed Keplerian circle 
at t = to- Blue dash-dotted line: osculating Keplerian circle after the first 
X axis crossing. As can be noted, its semimajor axis and eccentricity are 
equal to those of the initial unperturbed circle. Note also that after 2 yr the 
planet has not yet reached the x axis as it would have done in absence of 
mass loss. 
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(d) = for e — )• 0. The constancy of the osculating semimajor axis is not in 
contrast with the true trajectory, as clearly showed by Figure |H Again, the 
true orbital period is larger than the Keplerian one which, contrary to the 
eccentric case, remains fixed. Since 'D{E) = for e = and -7^(27r)|g^Q = 
—2iT, J^(0)|g^Q = 0, the radial shift per revolution is 

A.(2.)U = -^„(^). (45) 

Also in this case the secular increase of the radial distance is present, as 
qualitatively shown by Figure O Concerning At, after 2n it is 

Ar(27r) = —a f ; (46) 
n \nj 

also in this case, the orbital period is larger than the unperturbed one. 



3 The general relativist ic case 

The field equations of general relativity are non-linear, but in the slow- 
motion (/3 = v/c <^ 1) and weak-field [U/c^ <C 1) approximation they get 
linearized resembling to the linear equations of the Maxwellian electromag- 
netism; here v and U are the magnitudes of the typical velocities and the 
gravitational potential of the problem under consideration. This scenario 
is known as gravitoelectromagnetism [24^ I25j . In this case the space-time 
metric is given by 



ds^ = ( 1 - 2^ c^dt^ + -{H -drjclt- [1 + 2^] 6ijdx'dx\ (47) 



c' 

where, far from the source, the dominant contributions to the gravitoelec- 
tromagnetic potentials can be expressed as 

$=^, Jf=-^. (48) 

Here J is the proper angular momentum of the central body of mass M 
and r is so that r ^ GMjc? and r ^ J/{Mc); c is the speed of light in 
vacuum. For a non-stationary source the geodesic equations of motion yield 



®This fact has been quantitatively proven by the observation of convective binary stars 
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Ar (AU) 




Figure 5: Difference Ar{t) between the radial distances obtained from the 
solutions of two numerical integrations with MATHEMATICA of the equa- 
tions of motion over 3 yr with and without the initial conditions are 
the same. Just for illustrative purposes a mass loss rate of the order of 10~^ 
yr~^ has been adopted for the Sun; for the planet initial conditions corre- 
sponding to a = 1 AU, e = 0.0 have been chosen. The cumulative increase 
of the Sun-planet distance induced by the mass loss is apparent. 
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among other terms, — /3*(3 — f3 )^^q, i = 1, 2, 3 which, to order 0{c 
reduces to 

11 V GM f G m\ V 



For /t < such a perturbing acceleration is directed along the velocity of the 
test particle. Although of no practical interest, being of the order of W~^'^ 
m s~^ in the case of a typical Sun-planet system with M/M = —9 x 10^"'^'^ 
yr~^, I will explicitly work out the orbital effects of (|49p : the effects of the 
temporal variations of J have already been worked out elsewhere |27t I26| . 
Also in this case I will use the Gauss perturbative case. Since the radial and 
transverse components of the unperturbed velocity are 

naesinf 

Vr = -^J—, (50) 
Vr = (51) 

the radial and transverse components of (I49p . evaluated onto the unper- 
turbed Keplerian orbit, are 

_ 3fi nesmE 

c2 (l-ecos£;)2' ^^^^ 

_ 3/2 nVl-e2 



After lengthy calculations they yield 



da \ 6/i 



1 , (54) 



/ de\ 3fifi /r, \ f-, 1 



(2-e) . (55) 



\dt/ c'^a^e^ Vl^^ 



Contrary to the classical case, now both the osculating semimajor axis and 
the eccentricity increase for /i < 0. It turns out that the pericentre and the 
mean anomaly do not secularly precess. Also in this case the inclination and 
the node are not affected The qualitative features of the motion with the 
perturbation (j49p are depicted in Figure [6] in which the magnitude of the 
relativistic term has been greatly enhanced for illustrative purposes. 
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Figure 6: Black continuous line: true trajectory obtained by numerically 
integrating over 3 yr the equations of motion perturbed by ()49p . The planet 
starts from the perihelion on the x axis. Just for illustrative purposes, a 
factor —"iiijc? of the order of 5 x 10~^ AU yr~^ has been adopted for the Sun; 
for the planet initial conditions corresponding to a = 1 AU, e = 0.8 have 
been chosen. Red dashed line: unperturbed Keplerian ellipse at t = to = ^p- 
Blue dash-dotted line: osculating Keplerian ellipse after the first perihelion 
passage. As can be noted, its semimajor axis and eccentricity are larger 
than those of the initial unperturbed ellipse. 
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4 Discussion of other approaches and numerical 
calculations 



Here I will briefly review some of the results obtained by others by comparing 
with ours. 

Hadjidemetriou in Ref. [Hj uses a tangential perturbing acceleration 
proportional to the test particle's velocity v, 

(56) 

and a different perturbative approach by finding that, for a generic mass loss, 
the semimajor axis secularly increases and the eccentricity remains constant. 
In fact, with the approach followed here it would be possible to show that, 
to first order in — to), {a) = —{jj,/fi)a and (e) = and that the true 

orbit is expanding, although in a different way with respect to ^ as depicted 
by Figure El in which the magnitude of the mass- loss has been exaggerated 
for better showing its orbital effects. However, it must be noted that a term 
like (I56p is inadmissible in any relativistic theory of gravitation because it 
violates the Lorentz invariance. Indeed, this fact is explicitly shown for 
general relativity by Bini et al. in Ref. |26| where the full equations of 
motion of a test particle in a non-stationary gravitoelectromagnetic field are 
worked out (see, eq. (14) of Ref. [26j). In deriving them it is admitted that, 
in general, ^ = ^{t,r), but no gravitoelectric terms like ()56p occur. Instead, 
dl]) is compatible with eq. (14) of Ref. [26] . 

Schroder and Smith in Ref. jl], by assuming the conservation of the 
angular momentum, derive the orbital expansion by means of equations 
valid, instead, for orbits with constant radius only, i.e. v'^/r = ^{t)/r'^ 
and L = vr. Then, they assume that not only v but also r vary and put 
v{t) = -^f ^{t)/r, which is, instead, valid for circular orbits of constant radius 
only, into L = v{t)r{t) = vr getting fj,{t)r{t) = fj,r, where in my notation r 
and refers to the initial epoch to- With such an approach they obtain an 
expanded terrestrial orbit up to about 2 times larger than mine. 

Noerdlinger in Ref. [3| , following Jeans \T2\ and Kevorkian and Cole [18J , 
assumes for the variation of a quantity identified by him with the semimajor 
axis the following expression 

a{t)fi{t) = afi : (57) 

thus, his semimajor axis gets larger. Note that such an equation is the same 
obtained by Ref. By assuming a variation of /x linear in time (I57p would 
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Figure 7: Black continuous line: true trajectory obtained by numerically 
integrating with MATHEMATICA the perturbed equations of motion in 
Cartesian coordinates over 2 yr ; the disturbing acceleration (|56p has been 
used. The planet starts from the perihelion on the X axis. Just for illustra- 
tive purposes, a mass loss rate of the order of — 10~^ yr"^ has been adopted 
for the Sun; for the planet initial conditions corresponding to a = 1 AU, 
e = 0.8 have been chosen. Red dashed line: unperturbed Keplerian ellipse 
at i = to = ^p- Blue dash-dotted line: osculating Keplerian ellipse after the 
first perihelion passage. As can be noted, its semimajor axis is larger than 
that of the initial unperturbed ellipse, while the eccentricity remaines con- 
stant. Note also that after 2 yr the planet has not yet reached the perihelion 
as it would have done in absence of mass loss. 
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yield an increase of a according to 

d = - 0^ a > 0; (58) 

cfr. with my i^. As a consequence of the constancy of = fi{t)a{t)[l — 
e(t)^] and of ([57|) he obtains that the eccentricity remains constant, i.e. 

e = 0; (59) 

cfr. with my (jl2p . Moreover, another consequence of (j57p obtained by 
Noerdhnger in Ref. 0] is that the Keplerian period increases as 



pKep 



/x(t) 



(60) 



cfr. with my (|16p . Should the quantities dealt with by Noerdhnger are to be 
identified with the usual osculating Keplerian elements, his results would be 
incompatible with the real dynamics of a test particle in the field of a linearly 
mass-losing body, as I have shown. The quantity obtained by us which 
exhibits the closest resemblance with ()58p seems to be the secular variation 
of Ar(27r) for e = 0. Apart from matters of interpretation, the quantitative 
results are different. Indeed, I obtain for the Earth a secular variation of 
the semimajor axis of —2 x 10~^ m yr~^ and a shift in the radial position 
along the fixed line of the apsides of +1.3 x 10~^ m yr~^, while Noerdhnger 
in Ref. [1] gets a secular rate of his semimajor axis, identified with the 
Astronomical Unit, of about +1 x 10"^ m yr~^. Note that Noerdhnger uses 
for the Sun M/M = —9 x 10~^^ yr"^ as in the present work. 

Krasinski and Brumberg in Ref. [2] deal, among other things, with the 
problem of a mass-losing Sun in the framework of the observed secular in- 
crease of the Astronomical Unit for which, starting with an equation of 
motion like ([2]), they obtain an equation like (I58p . A mass- loss rate of 
—3 X 10~^^ yr~^j considered somewhat underestimated by Noerdhnger 
yields an increase of the Astronomical Unit of 3 x 10^'^ m yr^^. With such a 
value for I would obtain a decrease of the semimajor axis of —7 x 10~^ 
m yr~^ and an increase in r of +4 x 10~^ m yr~^. 

Concerning the observationally determined increase of the Astronomical 
Unit, more recent estimates from processing of huge planetary data sets by 
Pitjeva point towards a rate of the order of 10~^ m yr~^ \28\ I29| . It may 
be noted that my result for the secular variation of the terrestrial radial 
position on the line of the apsides would agree with such a figure by either 
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assuming a mass loss by the Sun of just —9 x 10^ yr" or a decrease of the 
Newtonian gravitational constant G/G ~ — 1 x 10^^'^ yr^^- Such a value for 
the temporal variation of G is in agreement with recent upper limits from 
Lunar Laser Ranging ^30j G/G = (2 it 7) x 10~^^ yr"^- This possibility is 
envisaged in Ref. [H] whose authors use a /a = —G/G by speaking about a 
small radial drift of —(6 it 13) x 10^^ m yr^^ in an orbit at 1 AU. I, now, 
apply my analysis to the daily decrement of the semimajor axis of LAGEOS 
[H] whose relevant orbital parameters are a = 12270 km, e = 0.0045. From 
([9]) it turns out that a secular decrease of a of the order of ~ 1 mm d^^ 
could only be induced by /ieZ/ie = —3 x 10^^ yi'^^- Clearly, it cannot be 
due to a variation of the Earth's mass M; if it had to be attributed to a 
variation of G, it would be orders of magnitude too large with respect to 
the bounds in Ref. [30] and Ref. [3T]. Adopting does not substantially 
alter the situation because the required /ie / ^0 would be only two orders of 
magnitude smaller. Now I look at the anomalous increase of the eccentricity 
of the lunar orbit amounting to CMoon = (0.9 it 0.3) x 10~^^ yr~^ [9]. Such a 
figure and ^ yield Ae//^e = 8.5 x lO"^^ y^-i. 

is one order of magnitude 
larger than the present-day bounds on G/G [30l [3T] . 

5 Conclusions 

I started in the framework of the two-body Newtonian dynamics by using a 
radial perturbing acceleration linear in time and straightforwardly treated 
it with the standard Gaussian scheme. I found that the osculating semima- 
jor axis a, the eccentricity e and the mean anomaly Ai secularly decrease 
while the argument of pericentre uj remains unchanged; the longitude of the 
ascending node and the inclination i are not affected by the phenomenon 
considered. The radial distance from the central body, taken on the fixed 
line of the apsides, experiences a secular increase Ar. For the Earth, such 
an effect amounts to about 1.3 cm yr~^. By numerically integrating the 
equations of motion in Cartesian coordinates I found that the real orbital 
path expands after every revolution, the line of the apsides does not change 
and the apsidal period is larger than the unperturbed Keplerian one. I 
have also clarified that such results are not in contrast with those analyt- 
ically obtained for the Keplerian orbital elements which, indeed, refer to 
the osculating ellipses approximating the true trajectory at each instant. I 
also computed the orbital effects of a secular variation of the Sun's mass in 
the framework of the general relativistic linearized gravitoelectromagnetism 
which predicts a perturbing gravitoelectric tangential force proportional to 
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v/r. I found that both the semimajor axis and the eccentricity secularly 
increase; the other Keplerian elements remain constant. Such effects are 
completely negligible in the present and future evolution of, e.g., the solar 
system. 

As a suggestion to other researchers, it would be very important to 
complement my analytical two-body calculation by performing simultaneous 
long-term numerical integrations of the equations of motion of all the major 
bodies of the solar system by including a mass-loss term in the dynamical 
force models as well to see if the N-body interactions in presence of such 
an effect may substantially change the picture outlined here. It would be 
important especially in the RGB phase in which the inner regions of the 
solar system should dramatically change. 
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